Eulerian and Lagrangian propagators for the adhesion model (Burgers dynamics) 
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Motivated by theoretical studies of gravitational clustering in the Universe, we compute propaga- 
tors (response functions) in the adhesion model. This model, which is able to reproduce the skeleton 
of the cosmic web and includes nonlinear effects in both Eulerian and Lagrangian frameworks, also 
corresponds to the Burgers equation of hydrodynamics. Focusing on the one-dimensional case with 
power-law initial conditions, we obtain exact results for Eulerian and Lagrangian propagators. We 
find that Eulerian propagators can be expressed in terms of the one-point velocity probability dis- 
tribution and show a strong decay at late times and high wavenumbers, interpreted as a "sweeping 
effect" but not a genuine damping of small-scale structures. By contrast, Lagrangian propagators 
can be written in terms of the shock mass function - which would correspond to the halo mass 
function in cosmology - and saturate to a constant value at late times. Moreover, they show a 
power-law dependence on scale or wavenumber which depends on the initial power-spectrum in- 
dex and is directly related to the low-mass tail of the shock mass function. These results strongly 
suggest that Lagrangian propagators are much more sensitive probes of nonlinear structures in the 
underlying density field and of relaxation processes than their Eulerian counterparts. 
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PACS numbers: 98.80.-k, 98.80.Bp, 98.65.-r, 47.27.Gs 



I. INTRODUCTION 

The formation of the large-scale structures observed 
in the present Universe is an important topic of modern 
cosmology In the standard scenario, these large-scale 
structures (such as galaxies, clusters of galaxies, filaments 
and large voids) have formed through the amplification 
by gravitational instability of small primordial Gaussian 
fluctuations, with a nearly scale-invariant power spec- 
trum, generated by an early inflationary stage. Within 
the concordant model of cosmology (see for instance [2| 
for comparison with current observational data), where 
the dynamics has been governed since z ~ 3000 by a 
collisionless dark matter component (although dark en- 
ergy dominates the expansion at about z < 1), this gives 
rise to a hierarchical evolution, as increasingly larger 
scales turn nonlinear (very large scales being in the linear 
Gaussian regime and small scales in the highly nonlinear 
regime). Today, scales beyond ~ 10 Mpc are well de- 
scribed by linear theory while scales below ~ 1 Mpc are 
within the highly nonlinear regime. 

The highly nonlinear regime has proved very difficult 
to handle by analytical tools so far (see for instance (H), 
and one must resort to numerical simulations and phe- 
nomenological models (such as the halo model [H) which 
involve some free parameters that are fitted to numer- 
ical results. However, weakly nonlinear scales remain 
within the reach of systematic analytical methods, based 
on various perturbative schemes. This regime is of great 
practical interest, as it is the focus of several observa- 
tional probes that aim at measuring the recent expan- 
sion history of the Universe and the growth of density 
fluctuations to constrain cosmology. This is the case for 
instance of baryon acoustic oscillation studies || 0] and 
weak lensing surveys Q. Moreover, on these large scales 



the system can be described through an hydrodynami- 
cal approach [![, which greatly simplifies the problem as 
compared with the full Vlasov equation that would be 
required to describe small scales where multi-streaming 
plays a key role [|[ . 

This has led to renewed interest in perturbation theory 
techniques, in order to increase the range where analyti- 
cal results apply. In the standard perturbative approach 
(see [3J for a review), the equations of motion for the 
density contrast and velocity fields, <5(x, t) and v(x, t), 
are solved as expansions over powers of the linear grow- 
ing modes <5l(x, t) and Vj v (x, t). Then, truncating these 
series at a given order and next performing the average 
over the Gaussian initial conditions yields the statistical 
quantities of interest, such as the density power spec- 
trum (which is the Fourier transform of the two-point 
density correlation) . As described in [E E3] , it is possi- 
ble to reorganize these perturbative expansions in terms 
of statistical quantities of interest, such as the two-point 
correlations and propagators, and to perform partial re- 
summations to improve the convergence. Alternatively, 
writing the equations of motion in terms of the two-point 
correlations and propagators themselves, one obtains an 
infinite hierarchy of equations (because one starts from 
nonlinear equations of motion) that must be truncated at 
some order. This can be done in many ways, using sys- 
tematic expansions of path-integral formulations of the 
dynamics (H El El or closure methods that truncate 
the hierarchy, written as a set of relations between the n 
and n + 1 correlations, by using some approximation at 
a level n + 1 to obtain a closed system of equations for 
lower-order quantities [3, EH • These methods can also 
be applied within a Lagrangian framework 

HE El, where 

the dynamics is written in terms of the displacement field. 
A simultaneous comparison of such approaches with nu- 
merical simulations is presented in [18j |. 
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Most of these alternative approaches to the standard 
perturbative expansion involve the so-called propagators, 
G(x, t;q), that can be seen as the cross-correlation be- 
tween a nonlinear field (f>(x, t) and its initial condition 
</>o(q) (or equivalently its linear growing mode <^>i(q, t)), 
as G(x, t;q) = (0(x, i)0o(q))- They are also restric- 
tions to £2 = of the more general response func- 
tions i2(xx, ti; X2, £2)1 that measure the mean change 
A0(xi, ti) of the field at time t\ generated by an infinites- 
imal perturbation A0(x 2 ,i2) at the earlier time £2 < t\- 
The partial resummation proposed in [jj [Io| predicts a 
Gaussian decays in the nonlinear regime, of the form 
e~ D + k where D + {t) is the linear density growth rate, 
which is in good agreement with numerical simulations, 
whereas some other approaches only yield a power-law 
decay [TTL [T^ . However, as explained in this strong 
damping is due to a "sweeping effect" rather than to a 
genuine loss of memory associated with relaxation pro- 
cesses. Thus, extending such resummation schemes to 
a Lagrangian framework [T3], one finds that Lagrangian 
propagators keep growing in the nonlinear regime instead 
of decaying (within these approximations). 

In order to shed light on the behavior of these response 
functions it is desirable to obtain their properties for 
some closely related dynamics where exact results can be 
derived. This would help understanding which physical 
processes govern their behavior and building more realis- 
tic approximations. The simplest dynamics that arises in 
this context is the Zel'dovich approximation (20j . where 
the particle displacement field is set equal to its linear 
prediction. Then, explicit expressions can be obtained 
for the correlation and response functions [^, and 
different-time Eulerian statistics show a Gaussian decay 
of the form e~ D + k at high wavenumbers, whereas La- 
grangian statistics remain equal by construction to their 
linear prediction. This simplified dynamics is useful to 
show how such a high-fc decay is produced. However, 
because its Lagrangian structure is so simple and it is 
grossly inaccurate in the highly nonlinear regime (par- 
ticles keep escaping to infinity after shell-crossing) it is 
insufficient to shed light on the strongly nonlinear regime 
and on the generic properties of both Eulerian and La- 
grangian statistics. 

A second dynamics, which is able to handle some non- 
linear effects in both Eulerian and Lagrangian frame- 
works, is provided by the adhesion model [2l|. This 
adds to the Zel'dovich dynamics an infinitesimal vis- 
cosity so that particles cannot cross. This binds col- 
lapsed structures and fixes (at least at a qualitative level) 
the main failure of the Zel'dovich approximation (but 
the dynamics of collapsed halos remains non-trivial, see 
|22l]). Numerical simulations show that this simple dy- 
namics already provides a significant improvement over 
the Zel'dovich dynamics and is able to reproduce the 
large-scale skeleton of the cosmic web [23| . On the other 
hand, as seen in [2l| , the adhesion model also corresponds 
to the standard Burgers equation 24J ) . which was origi- 
nally introduced as a simplified model for hydrodynami- 



cal turbulence and has been the subject of many studies, 
see [25j for a recent review. The advantage of this nonlin- 
ear dynamics is that it can be explicitly integrated, and 
exact results can be derived in some cases, mostly in one 
dimension for power-law initial conditions. 

Therefore, in this article we investigate the proper- 
ties of the Eulerian and Lagrangian propagators, or re- 
sponse functions, R(x,t;qo), obtained within this adhe- 
sion model (we do not consider here equal-time power 
spectra or correlation functions, which have been stud- 
ied in previous works, e.g. [26|, H3]). Since we are in- 
terested in exact analytical results, we focus on the one- 
dimensional case for power-law initial conditions, where 
asymptotic tails can be derived in both linear and highly 
nonlinear regimes. We also consider in more details the 
two representative cases of Brownian and white-noise ini- 
tial velocity, where explicit expressions can be obtained 
for the full propagators. We first recall in section |TT] how 
the Zel'dovich and Burgers dynamics can be derived from 
the cosmological gravitational dynamics and we present 
the one-dimensional power-law initial conditions that we 
consider in this paper. We also recall how the Burgers 
dynamics can be integrated through the Hopf-Cole trans- 
formation [28, 29[ and its geometrical interpretation, and 
we describe the self-similar evolution that is obtained for 
these power-law initial conditions. Then, we study the 
Eulerian propagators in section IIII1 and we recover the 
"sweeping effect" described above. Next, we consider the 
Lagrangian propagators in section IIVI We derive their 
relation with the density field (through the shock mass 
function) and find as expected that their properties are 
quite different from their Eulerian counterparts. Finally, 
we conclude in section fVl 



II. INITIAL CONDITIONS AND 
GEOMETRICAL SOLUTION 

A. The adhesion model 

We briefly recall here how the Burgers equation ap- 
pears in the cosmological context, through the "adhesion 
model" [U |3(J. On scales much larger than the Jeans 
length, both the cold dark matter and the baryons can 
be described as a pressureless dust. Then, neglecting or- 
bit crossings one can use a hydrodynamical description 
governed by the equations of motion [l| , 



g + V.[(l + £)v] 
— + Hv+ (v.V)v 

OT 



= 0, 



A<J) = 



-V4>, 



(i) 
(2) 

(3) 



where r = J dt/a is the conformal time (and a the scale 
factor), Ti. = dlna/dr the conformal expansion rate, and 
fl m the matter density cosmological parameter. Here, S 
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is the matter density contrast and v the peculiar velocity. 
Since the vorticity field decays within linear theory [l|, 
the velocity is usually taken to be a potential field, so 
that v is fully specified by its divergence, —9, or by its 
potential, ip, with 



= -V.v, 



-V?/>, whence 6 = Aip. (4) 



In the linear regime, one finds that the linear growing 
mode satisfies (using a subscript L for linearized quanti- 
ties) 



9l = fH5h whence c/)l — 



IpL, 



(5) 



where /(r) is defined from the linear growing rate D + (t) 
of the density contrast by f — dlnD + /dlna, and D+(t) 
is the growing solution of 



(6) 



If we make the approximation that relation ([5]) remains 
valid in the nonlinear regime, that is, we replace the Pois- 
son equation ([3]) by the second Eq.([5]), <j> — 3f2 m W"0/(2/), 
then we obtain for the Euler equation 



dv 

7h 



— + i 



0. 



(7) 



Obviously, as shown by Eq.0, within this approxima- 
tion the velocity field now evolves independently of the 
density field. As is well known plj . approximation 0) is 
actually identical to the Zel'dovich approximation [20(. 
Indeed, a change of variables for the velocity field yields 



9u , _, n ( dDj 

+ (u.V)u = with v = ' 



3D 



\ At 



(8) 



Equation ([5]) is the equation of motion of free particles, 
du/dl? + = 0, hence the trajectories are given by 



x = q + £» + (r)u L0 (q), v 



ULoiqj 



(9) 



where q is the Lagrangian coordinate and s = D+u^o 
is the displacement field that is exactly given by the lin- 
ear theory. Equation ([9]) is the usual definition of the 
Zel'dovich approximation (i.e. setting s = Sz,). Then, in 
order to prevent particles from escaping to infinity after 
shell-crossing, and to mimic the trapping within grav- 
itational potential wells, one can add an infinitesimal 
viscosity, vAu, to the right-hand-side of Eq.©, which 
becomes the standard Burgers equation [24[ • This is the 
adhesion model proposed in (2lj | to study the formation 
of large-scale structures. 

Numerical simulations [23| have shown that this sim- 
plified model is able to reproduce the cosmic web seen in 
gravitational simulations: starting with identical initial 
conditions it recovers the shape and the location of fila- 
ments (but the latter are now infinitesimally thin while 



halos are point-like objects). From a theoretical point 
of view, the advantage of the adhesion model and of the 
Zel'dovich approximation is to provide a simpler dynam- 
ics which can be exactly solved (at least in a few non- 
trivial cases) while remaining close to the gravitational 
dynamics. In particular, the advective quadratic nonlin- 
earities that appear in the equations of motion (jX|)- (J3J) are 
preserved. Then, both Zel'dovich and Burgers dynam- 
ics may be used as benchmarks to test approximation 
schemes devised for the gravitational dynamics [l^, [3l[ . 

For the Zel'dovich dynamics it is possible to obtain 
explicit expressions for many quantities of cosmological 
interest, such as the nonlinear matter power spectrum 
[32l HU and the propagators [t| Ojj|. They all show a 
Gaussian damping at high wavenumbers due to the un- 
bounded random displacement of particles which erases 
all structures in the highly nonlinear regime (as particles 
keep traveling with their initial velocity u after shell- 
crossing). This motivates the study of the Burgers dy- 
namics, which avoids this spurious damping. Moreover, 
while the Zel'dovich dynamics is ill-defined for initial con- 
ditions with a slope — 1 < n < 1 (because of the strong 
power at large wavenumbers the linear displacement field 
shows UV divergences), both the Burgers and the grav- 
itational dynamics remain well-defined (as viscous stick- 
ing or gravitational trapping regularize the dynamics on 
small nonlinear scales). Since the range of cosmological 
interest is — 3 < n < 1 it is useful to study a model that 
covers this domain. 

However, whereas both Zel'dovich and Burgers dynam- 
ics can be explicitly integrated, the expression obtained 
for the latter is much more complex, see Eqs.p5 |) - (p?i?)) 
below. As applications in cosmological context require 
the computation of ensemble averages over the initial 
conditions, obtaining explicit results when the relation 
between the initial conditions and the final result is in- 
volved is all the more difficult. It turns out that it is only 
possible to derive exact results for a few specific cases, in 
one dimension. This is the reason why in the following 
we consider the one-dimensional Burgers equation, (jXT 
and focus on power-law initial conditions. 



B. Equations of motion and initial conditions 

As explained above, we are led from Eq.© to the study 
of the one-dimensional Burgers equation. Thus, making 
the change of notation D+ — ► t for simplicity, we obtain 
the standard one-dimensional Burgers equation for the 
velocity field u(x, t) in the limit of zero viscosity, 



du 
dt 



du 
dx 



d 2 u 
dx 2 



with 



H 



(10) 



while the density field still obeys the usual continuity 
equation (see [22| for more details), 



dp 
dt 



with 



p(x,0)=p o . (11) 
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Thus, the initial conditions are set at t = (which cor- 
responds to D + = 0), with a uniform density p (which 
corresponds to the mean comoving matter density) and a 
Gaussian random velocity field uo(q). Introducing again 
the velocity divergence and potential as 

dx ' dx dx 2 ' 

and normalizing Fourier transforms as 

/oo 
dke ikx 6(k,t), (13) 
-oo 

the initial divergence, —9q, is taken as Gaussian, homo- 
geneous, and isotropic, so that it is fully described by its 
power spectrum Pg (k) with 

{Oo) = 0, (0o(*i)0o(*2)) = Soikt+^PgM, (14) 

where So is the Dirac distribution. As in [3lj . but re- 
stricting ourselves to one dimension, we consider power- 
law initial power spectra, 

P 8o (k) = — k n+2 with - 3 < n< 1, (15) 

where D is a normalization factor. Since we have 
u(k,t) — {i/k)6(k,t), the initial energy spectrum is a 
power law, 

(fio(fci)tio(fe)) = M*i + k 2 )E (kx), (16) 

with 

E (k)=k- 2 Pe (k) = ^-k n . (17) 

As can be seen from the analysis in .31], the in- 
dex n introduced in Eqs. (fll^ - (fr7|) also corresponds to 
the standard index n used in three-dimensional cosmol- 
ogy, where it is defined from the linear density contrast 
Sl as P^~ 3 (fc) cx k n . More precisely, in arbitrary di- 
mension d one must define the initial power spectra as 
Pe (k) cx k n+3 ~ d and E (k) cx k n+1 - d . Then, many 
properties only depend on the index n, independently of 
dimension d, such as the scaling laws ((20)) and (|38|) below, 
see [13 . 

Since for the standard CDM cosmology the local slope 
n runs from 1 at large scales to —3 at small scales, the 
range (fT5|) covers the cases of cosmological interest. This 
interval can actually be split into two distinct classes. 
First, for — 1 < n < 1, which corresponds to large power 
at high k ( "UV-class" ) , the initial velocity field is ho- 
mogeneous. Moreover, it is singular (e.g., a white noise 
for n — 0) but this ultraviolet divergence is regularized as 
soon as t > by the infinitesimal viscosity [24J . Note that 
this is a non-perturbative effect and that the Zel'dovich 
dynamics, which lacks this regularization process, is not 
defined in this case. This class is the "type B initial con- 
ditions" studied in 34]: the initial velocity is the deriva- 
tive of a fractional Brownian motion. 



Second, for —3 < n < —1, which corresponds to large 
power at low k ("IR-class"), the initial velocity field is 
no longer homogeneous but only shows homogeneous in- 
crements [35[ (thus the divergence, —6, is still homoge- 
neous) . This is the "type A initial conditions" of [34[ : the 
initial velocity is a fractional Brownian motion. Then, 
to handle the infrared divergence at low k one must 
choose a reference point, such as the origin xq = 0, with 
uo{xq) = 0, and define the initial velocity in real space 
as 

/oo 
dfc (e ikx - e ikxo ) u (fc), for -3 < n< -1. 
-OO 

. (18) 

Note that because of the nonlinear advective term in the 
Burgers equation ([TO]) , the increments of the velocity field 
are no longer homogeneous for t > 0, which also means 
that the divergence 9(x,t) is no longer homogeneous ei- 
ther. However, at large distance from the reference point 
(i.e. taking the limit \xq\ — > oo or |x| — > oo), we can 
expect to recover an homogeneous system (in terms of 
velocity increments and matter distribution), see [36j for 
more detailed discussions. This can be shown explicitly 
for the case n = —2, where the initial velocity field is a 
Brownian motion [2y, [3?| • On the other hand, we may 
add a low-fc cutoff to the initial power spectrum and re- 
strict ourselves to finite times and scales where the in- 
fluence of the infrared cutoff is expected to vanish for 
equal-time statistics. This property can also be seen at 
a perturbative level or from numerical simulations in the 
three-dimensional gravitational case studied in cosmol- 
ogy for equal-time statistical quantities (38l . |39| . This is 
due to the Galilean invariance of the equations of motion 
and also holds for the Zel'dovich dynamics [l9j . 

Normalizing the velocity potential by ipo(Q) = 0, that 

is, 

M%)= I dx'u (x'), (19) 
Jo 

and the velocity field by u (0) = if — 3 < n < —1, the 
initial conditions obey the scaling laws 

A>0 : 6 (\x) ^ \-( n+ ^ 2 o (x), 
u (Xx) ^ A-(" +1 )/ 2 u (x), 
MAx) '= A (1 -" )/2 Vd(aO, (20) 

where " = " means that both sides have the same statis- 
tical properties. 

These initial conditions can also be expressed in terms 
of the linear density field as follows. If we linearize the 
equations of motion (P~CT|) - (| 1 1 [) we obtain the solution for 
v = 0, 

B L (x,t)=0 o {x), 5 L (x,t)=te (x), (21) 

where we defined the density contrast as 5(x, t) = 
(p(x,t) — po)/po and the subscript L stands for the linear 
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regime. (Note that we recover the usual growing mode for 
the density contrast since as explained above t stands for 
the linear growing mode D + in the cosmological context.) 
Therefore, we can as well define the initial conditions by 
the linear density contrast <5r,(x,i), which is Gaussian, 
homogeneous, and isotropic, with a power spectrum 



3<n<l: Ps L (k,t) =r Pe (k) oar k 



2 7,1+2 



(22) 



This is the manner in which initial conditions are usually 
defined in the cosmological context. (As explained above, 
in d dimensions we would obtain Pg L (k,t) oc t 2 k n + 3 ~ d 
which recovers the cosmological notation for d = 3.) 

In this article we consider in more details two repre- 
sentative cases where many exact results can be obtained 
[Hii^lli, [13, SO, EH SI. First, the case rc = 0, asso- 
ciated with the UV-dominated range —1 < n < 1, corre- 
sponds to a white-noise initial velocity field, normalized 



n = 0: (u o (q))=0, (u (qi)u (q 2 )) = D 8 D {qi - q 2 ), 

(23) 

(V>o(?))=0, (tpo(qi)il>o{q2)) = Dqi, for < q x < q 2 . 

(24) 

Thus, the initial velocity potential is a bilateral Brownian 
motion that starts from the origin. 

Second, the case n = —2, associated with the IR- 
dominated range —3 < n < — 1, corresponds to a Brow- 
nian initial velocity field, normalized as 

n = -2: (e o (q))=0, (6 (qi)0 (q2)) = D8 D {q x - q 2 ), 

(25) 

(wo(<?))=0, (u (qi)u (q 2 )) = D qi, for < qy < q 2 . 

Thus, it is now the initial velocity field which is a bilateral 
Brownian motion that starts from the origin. 



C. Hopf-Cole solution and first-contact parabolas 



As is well known [28, 29], making the change of variable 
ip(x,t) = 2fhx'E(x,t) transforms the nonlinear Burgers 
equation into the linear heat equation. This gives the 
explicit Hopf-Cole solution 



ip(x,t) = 2vln 



dq 



exp 



ip (q) {x-qf 



Then, in the inviscid limit v 
method gives 



2v Avt 

(27) 

+ the steepest-descent 



i/>(x, t) — max 



(x - qf 



2t 



,(28) 



u(x,t) = u (q) = X q ^ t \ (29) 

where we introduced the Lagrangian coordinate q(x, t) 
defined as the point where the maximum in Eq. ()28|) is 



reached. In particular, this is the Lagrangian coordi- 
nate (i.e. the initial location) of the particle that is lo- 
cated at the Eulerian position x at time t. Here and 
in the following we note by the letter q the Lagrangian 
coordinates, which appear in the arguments of the ini- 
tial fields at t — 0, and by the letter x the Eulerian 
coordinates, which appear in the arguments of the Eu- 
lerian fields at any time t > 0. The Eulerian locations 
x where there are two solutions, g_ < to the maxi- 
mization problem ([2"5]) correspond to shocks and all the 
matter initially between q_ and q + is gathered at a;. At 
these points the velocity is discontinuous while the den- 
sity is infinite. The application q i— > x(q,t) is usually 
called the Lagrangian map, and x i— > q(x, t) the inverse 
Lagrangian map (which is discontinuous at shock loca- 
tions) [25j. Outside of shocks Eq. (|2"9")) clearly shows that 
one recovers the free-streaming dynamics ((5J -([9]) in the 
inviscid limit v —> + . 

The maximization problem (|28p has a well-known geo- 
metrical solution [24[. Indeed, let us consider the family 
of upward parabolas P x ,c{q) centered at x and of height 
c, with a curvature radius t, 



v x ,M = 



(g - xf 
2t 



(30) 



Then, moving down V XtC (q) from c = +oo, where the 
parabola is everywhere well above the initial potential 
V'o(g) [ill, until it touches the curve ipo(q), the abscissa 
q of this first-contact point is the Lagrangian coordinate 
q(x,t). If first-contact occurs simultaneously at several 
points there is a shock at the Eulerian location x. One 
can build in this manner the inverse Lagrangian map 

x i — ► q(x, t). 

Finally, the continuity equation (fTTj) can also be inte- 
grated as follows in one dimension (see [22j for a discus- 
sion of the more complex case of higher dimension) . The 
conservation of matter implies that the density field is 
related to the inverse Lagrangian map, x i— > q(x,t), and 
to the velocity potential, ip(x,t), through the Jacobian 



p(x,t)dx = p dg, 



(31) 



whence, using Eq. 



p(x, t) = po — = p 



1 -t- 



du 
dx 



Po 



l + t 



dx 2 



which gives for the matter density contrast, 

d 2 ib 

s(x,t) =te[x,t) =t-g^- 



(32) 



(33) 



Here we used the fact that particles do not cross each 
other, so that x(q) and q(x) are monotonous increasing 
functions and there is no need to keep the absolute value 
for the Jacobian J = \dq/dx\. In the last two equal- 
ities in (|3l?|) we used Eg. (f2"9")l and the definition of the 
velocity potential. One can easily check that ([3"!?]) is also 
valid for shocks, which give rise to Dirac density peaks. 
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We can note from the last expression (|33[) and Poisson's 
equation that the velocity potential, ijj, is equal to the 
gravitational potential, <fi, up to a normalization (and an 
additive quadratic term oc pox 2 /2 associated with the 
mean density po). This property, which only holds at the 
linear order in higher dimensions, is also at the basis of 
the Zel'dovich approximation used in cosmology and re- 
called in section HTX1 sec Eq.© (see also the discussion 
in section 2.2.2 of [301). Thus, Eq. ([33"|) recovers the well- 
known fact that the Zel'dovich approximation is actually 
exact in one dimension before shell-crossing. 

For the UV-class, — 1 < n < 1, shocks are expected to 
be isolated and in finite number per unit length, whereas 
between shocks the inverse Lagrangian map, x t— > q(x, t), 
is constant, so that the density field is made of a finite 
number of Dirac peaks per unit length amid empty space. 
For the IR-class, —3 < n < — 1, shocks are expected to be 
dense so that the density field is made of an infinite num- 
ber of Dirac peaks per unit length (without any smooth 
background). These results have been explicitly proved 
for white-noise [U |43[ and Brownian [42j initial velocity; 
they are only supported by phcnomenological arguments 
and numerical simulations [34J for generic values of n. 



D. Self-similar evolution 



For the initial conditions (|15[) . using the scaling laws 
(|2U|) one can see from the explicit solution (|28p that the 
nonlinear Eulcrian fields obey the scaling laws 

4>(X, t) 'l= W id-«)/(n+3) ^ (t-2/(n+3) Xj ^ > (34) 

u(x,t) = w t-(«+i)/(»+3) u (V 2 /("+ 3 ).T,l) , (35) 
q(x, t) ^ t 2 /("+ 3 ) q (t~ 2 ^ n +^x, l) , (36) 
5(x,t) '= 8{t- 2 '^x,l) . (37) 



We can check that the scaling ([37]) agrees with the linear 
mode ([2T]) . These scalings mean that the dynamics is 
self-similar: a rescaling of time is statistically equivalent 
to a rescaling of distances, as 



A > : t 



Xt. 



\ 2 '^x. 



(38) 



Thus, as in the standard cosmological scenario the 
system displays a hierarchical evolution as increasingly 
larger scales turn nonlinear. More precisely, since in the 
inviscid limit there is no preferred scale for the power-law 
initial conditions (fTS"]) , the only characteristic scale at a 
given time t is the so-called integral scale of turbulence, 
L(t), which is generated by the Burgers dynamics and 
grows with time as in (|38p , 



(39) 



L(t)<xt 2 ^ n+3 \ 



It measures the typical distance between shocks, and it 
separates the large-scale quasi-linear regime, where the 



energy spectrum and the density power spectrum keep 
their initial power-law forms, from the small-scale non- 
linear regime, which is governed by shocks, where the 
density power spectrum reaches a universal white-noise 
behavior (i.e. Pg(k,t) has a finite limit for k ^> 1/L(t)) 
[UHSHl. Note that the scalings (|33] ) -(|39" |) hold for any 
dimension d, provided we define the initial conditions by 
Ps L {k,t) cx t 2 k n+3 - d . 

In order to express the scaling laws (j34"|) - (j37| it is con- 
venient to introduce suitable dimensionless scaling vari- 
ables. 



Q 



L(t) 



x tu 



L(t) 



L(ty 



(40) 



where L(t) is the characteristic scale ([31?]). which we nor- 
malize as 



L(t) = (2DT) 



2Nl/(n+3) 



(41) 



where the constant D was introduced in Eq. (|15p . Thus, 
equal-time probability distributions written in terms of 
these variables no longer depend on time, and the scale 
X = 1 is the characteristic length of the system, at any 
time. On large quasi- linear scales, I>1, density fluctu- 
ations are small and the distributions are strongly peaked 
around their mean, with tails that are directly governed 
by the initial conditions (but shocks cannot be neglected 
for n > —2). On small nonlinear scales, I < 1, den- 
sity fluctuations are large and probability distributions 
show broad power-law regions [26j, |27J, [3l| . Note that for 
the associated power-law initial conditions in the three- 
dimensional cosmological context, the gravitational dy- 
namics also develops the same self-similar evolution (|39[) , 
see [l|. 



III. EULERIAN PROPAGATORS 

A. Relation with the velocity probability 
distribution 

We now consider the Eulcrian propagator of the ve- 
locity potential, R^(x,t;qo), defined as the functional 
derivative of ij;(x,t) with respect to the initial potential 
^0(90) at point q , 



R v {x,t;q ) = ( — — ; 

Vip (qo) 



(42) 



Then using the explicit Hopf-Cole solution l|27p to per- 
form the functional derivative and after the inviscid limit, 
v — > + , has been taken, we have 



R*(x,t\qo) = (5 D [q(x,t)-q }). 



(43) 



Here q(x, t) is again the inverse Lagrangian map intro- 
duced in (|29|) . From its functional definition we can see 
that the Eulerian propagator RV(x,t;qo) describes the 
sensitivity of the nonlinear potential at a given time t 
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with respect to the initial conditions. This response func- 
tion can also be seen as a memory kernel, with a time- 
dependence that would give an estimate of the time-scale 
beyond which initial fluctuations at a given wavelength 
appear to be damped. Note that from the expression (|43|) 
it obeys the sum rule 



dq R*(x,t;q ) = 1. 

Next, the average (|42|) gives 

R?(x,t;q ) =p x (qo,t), 



(44) 



(45) 



where p x (q, t) is the one-point probability distribution of 
the Lagrangian coordinate q(x,t). Since shocks form a 
set of zero measure, we can use (f2l?f to write Eq. (|4"5|) as 



R?(x,t;q ) = -p x (u, t), with 



x - go 
t 



(46) 



where p x (u,t) is the one-point Eulcrian velocity proba- 
bility distribution. Then, for the IR-class —3 < n < — 1 
where the velocity field is not homogeneous, as discussed 
in section ITTBl the propagator R^(x, t;q ) is not homoge- 
neous either, whereas it is homogeneous for the UV-class 
— 1 < n < 1. In the latter case, it can be useful to go to 
Fourier space, with the normalization 

-Kn<l: (H^-) = S D (k-k )R^(k,t), (47) 

where the Dirac factor Spik — ho) is due to the statistical 
invariance through translations. This yields 



l<n<l: R^(k,t) 



dxe~ ikx R^{x,t;0). 



Within the cosmological context, where one is mostly 
interested in the density field, one rather considers the 
density propagator R s defined by 



R?(x,t;q ) = ( 



V5{x,t) , 
T>S L0 (q )' 



(49) 



where 6(x, t) = {p(x, t) — po)/po is again the density con- 
trast and from Eq. lf2"Tj) we defined <$lo(<7o) as 

Slo (qo ) = #o0?o) whence S L (x,t) = tS L0 (x). (50) 

Then, from Eqs. ((50|) and (l33|) . we have in Fourier space 

S LQ (k) = -k 2 %(k) and S(k,t) = -tk 2 $(k,t). (51) 

Therefore, the density propagator (|4T)|) is equal to the ve- 
locity potential propagator defined in Eq. (|42p , multiplied 
by time t, 

R s (k.t) = tR^(k,t) and R s (x, t; q ) = tR*(x, t; q ). 

(52) 



For the power-law initial conditions ([15]) the first Fourier- 
space equality in (|52[) only holds for — 1 < n < 1 as in 

It is interesting to note that Eqs. ([46|) and (|52| show 
that Eulerian propagators, or response functions, are un- 
likely to be effectual in probing the nonlinear structures 
built in the underlying density field. As discussed in 
the following sections, the relation (|46f implies a strong 
dependence on a "sweeping effect" associated with long 
wavelength modes of the velocity field. Besides, whatever 
the magnitude of this effect, it is clear that the one-point 
velocity distribution does not provide much pertinent in- 
formation on the density field. Therefore, in more general 
dynamics such as gravitational clustering, which should 
still behave in a similar fashion, computing the Eulerian 
propagators in a perturbative manner up to high order, 
or exactly in cases such as the Burgers dynamics stud- 
ied here, is unlikely to shed much light on the density 
field. This is a strong motivation to study Lagrangian 
propagators, which are much more directly linked to the 
properties of the underlying density field as we shall find 
out in Sect. IIVI below. 



B. Linear regime and IR-class 

At large scales and early times density and velocity 
fluctuations are small and the system is well described 
by the linearized equations of motion. Note however that 
such a linear regime does not exist for the UV-class — 1 < 
n < 1, where the initial velocity variance (u 2 ) shows a UV 
divergence and shocks dominate the dynamics as soon as 
t > 0. For —2 < n < — 1 shocks also play a key role as 
soon as t > 0, and modify the naive linear predictions at a 
quantitative level (i.e. numerical prefactors), but they do 
not change the qualitative behavior of the distributions 
seen at large scales or early times (i.e. exponents in the 
exponential tails), see (3l| . 

Thus, for the IR-class, or for generic initial conditions 
with a high-fc cutoff, it is interesting to consider the linear 
predictions for the Eulerian propagators (|42|) and (|49|) . 
From Eq. (|2ip we have ipL{x, t) = ipo{x) which would give 
at zeroth-order R^(x,t;qo) = 5d{x — qo). However, we 
can obtain the first-order prediction from Eq. (|46p . In- 
deed, for Gaussian initial conditions the initial velocity 
probability distribution at position q reads as 



P q (uo) 



1 



2tt<7 Uo (q) 



,-"o/(2<(9)) 



where we introduced the initial velocity variance 
<(q) = {uo(qf). 



(53) 



(54) 



For initial conditions with an IR cutoff the initial velocity 
field can be homogeneous, so that a 2 Q and p q (u ) do not 
depend on position q, but for the IR-class —3 < n < — 1 
this is not the case, as discussed in section IIIBI above 
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Eq. (fT8|) . Then, the linear prediction for the Eulerian 
propagators reads as 



Rt(x,t; qo ) = _^_ e -(— )"/(«"<(•)), (55) 
V2Trta Uo (x) 



R 5 L {x,t;q ) 



2na Uo (x) 



e -(x-« ) 2 /(2i 2 a^ (x))_ /ggs 



Of course, in the limit t — > we recover R^(x, 0; go) = 
(5_d(x - q ) and i? 5 (x,0;g o ) = 0. 

If the initial conditions show an IR cutoff, so that the 
system is homogeneous, we can go to Fourier space as in 
Eq.gZ]). This yields 



i^(M)=e- t2fe2 < /2 , R 5 L (k,t)=te~ t2k2 </ 2 . (57) 

We can note that for the Zel'dovich dynamics particles 
exactly follow the linear displacement field as recalled 
below Eq.©. Then, Eqs. (|55|) - ([57)l are exact (for homo- 
geneous systems) || [lj| and Eulerian propagators show 
the characteristic Gaussian decay e _< k due to the ran- 
dom advection of particles that eventually erases small- 
scale structures as particles escape to infinity after shell- 
crossing. For the adhesion model, that is the Burgers 
dynamics studied in this article, Eqs. (|55)) - (|57p only hold 
in the linear regime. 

If the initial conditions belong to the IR-class, —3 < 
n < — 1, without IR cutoff, the initial velocity variance 
depends on position q as 



3 < n < -I 



(q) = D> \q\-t n+1 \ (58) 



as seen from Eq. (j2"0)) , with D' = D for n = — 2. Then, 
at any finite time t, as we go far from the reference point 
xq = the nonlinear velocity distribution becomes dom- 
inated by the initial Gaussian distribution (as nonlin- 
ear effects have only redistributed matter over the finite 
scales X ~ 1 and U ~ 1), 



— 3 < n < —2, |x| — > oo 

Px(u) 



|(n+l)/2 



-u 2 \x\* l+1 /(2D') 



(59) 



For — 2 < n < — 1 (in one dimension) shocks modify the 
numerical factor in the exponential decay, but not the 
exponents that remain of the form hip x (u) ~ — |x| n+1 w 2 
[3j . Note that the asymptotic behavior (fB"9"| can be ex- 
plicitly derived from the exact velocity probability dis- 
tribution for the Brownian case n — — 2 [26]. This can 
also be related to the "principle of permanence of large 
eddies" in the generic case [46] . 

Then, from expression (|46[) we can see that the Eu- 
lerian propagators vanish as soon as t > in the limit 
|x| — » oo, that is far from the reference point, at fixed 
separation |x — go I- More precisely, in agreement with 
the sum rule (|44[) , the Eulerian propagator spreads over 
an increasingly large region so that it vanishes for any 
fixed distance |x — go I- Indeed, as we go far from the 



reference point to avoid boundary effects, the local ve- 
locity variance becomes increasingly large, because of 
long-wavelength modes, and this large collective veloc- 
ity leads to a very fast damping of Eulerian propaga- 
tors (that becomes instantaneous in the limit |x| — > oo) 
as structures are transported over large distances. This 
"sweeping effect" , which agrees with the previous dis- 
cussion below Eq. (|57p . is not a true loss of memory if 
the system is observed on a global scale, as it does not 
imply that small-scale structures are erased (thanks to 
the Galilean invariance of the equations of motion) since 
this divergence arises from low-fc modes and is associated 
with almost uniform random translations. 

This motivates the study of Lagrangian propagators, 
introduced in [I7| in the cosmological context, and con- 
sidered in section IIVI below, to go beyond this sweeping 
effect. We can note that the sensitivity of Eulerian prop- 
agators to long-wavelength modes of the velocity field 
can be directly seen from Eq. (|46[) . which is not invariant 
through Galilean transformations. 



C. White-noise initial velocity 

Before turning to Lagrangian propagators, we consider 
the case n = of white-noise initial velocity, that is rep- 
resentative of the UV-class. For such initial conditions, 
— 1 < n < 1, the initial variance a^ is infinite as it shows 
a UV divergence and the initial velocity field u (q) is not 
a regular function (unless we add a UV cutoff). Then, 
shocks govern the dynamics as soon as t > and the 
linear regime described in Eqs.([55 ]l -(|57 | does not exist 
[HJ|. Thus, it is difficult to derive generic results and we 
focus below on the white-noise case n = where exact re- 
sults can be obtained. This should provide a qualitative 
illustration of the behaviors obtained for — 1 < n < 1. 

In terms of the dimensionless scaling variables (|4H|) the 
Eulerian propagator (|45|) writes as 



R^(x,t; qo ) = ^-R^(X 



■Qo), 



with 



R*(X - Qo) = P(U) and U = X - Q , 



(60) 



(61) 



where P{U) is the one-point probability distribution of 
the dimensionless velocity U. The latter was obtained in 
S3 as 



P(U) = J(U)J(-U) 



-zoo 



ds 



2tt\ Ai(s) ' 
(62) 

where Ai(x) is the Airy function. Of course, because 
of the statistical homogeneity and isotropy of the sys- 
tem, we can check that R^(x, t; q ) only depends on the 
distance |x — go| and on time, through the combination 
\X — Qo\ thanks to the scale-invariance (I3"51) . Then, 
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log(|X Q |) log(K) 



FIG. 1: Upper panel: The Eulerian propagators R^(X — 
Qo) = R S (X — Qo) obtained for n = (white-noise ini- 
tial velocity), in terms of dimensionless variables. They are 
also equal to the one-point velocity distribution P(U), with 
U — X — Qo, from Eq. (|6ip . The dashed lines are the asymp- 
totic cubic exponential behavior (|63[) . Lower panel: Same as 
upper panel but on a logarithmic scale. 



FIG. 2: Upper panel: The Eulerian propagators R^(K) = 
R 5 (K) in Fourier space, obtained for n = (white-noise initial 
velocity), in terms of the dimensionless wavenumber K, from 
Eq. (|65p . The dashed line is the asymptotic 3/2-exponential 
behavior (I67[l for the amplitude of R/ , that is setting the 
cosine to unity. Lower panel: Same as upper panel but on a 
logarithmic scale. The cusps correspond to changes of sign 

for py. 



Eq.([S2]) gives the asymptotic behavior at large dimen- 
sionless separation X [23, 



\X\ > 1 : R*(X) 



21X1 



Ai'(-wi) 



,-^\X\-\X\ 3 /3 



(63) 



where —u}\ is the first zero of the Airy function [u)\ ~ 
2.338). We display in Fig. [1] the Eulerian propagator 
R^(X — Qo) = P(U) in terms of dimensionless variables. 
Going to Fourier space, Eqs.(f48|) and (f62|) yield 

^{k, t) = I&(K), with K = L(t) k = (2Dt 2 ) 1/3 k, 

(64) 

where we introduced the dimensionless scaling wavenum- 
ber K, and 



-ZOO 



1 



2th Ai(s + Lff/2)Ai(s - iK/2) ' 



(65) 



We can check that at low wavenumber R^(0) = 1. 
Indeed, on large scales particles have not significantly 
moved from their initial positions (at a finite time t par- 
ticles have moved over distances of order (2Dt 2 ) 1 ^ 3 , i.e. 
X ~ 1) so that R^(x,t;qo) ~ 6d(x — q) as seen from 



large distances or in the limit t — ► 0. More precisely, 
R^(x, t; qo) decays faster than any power-law over \x— qo\ 
at large distances, so that R^(k, t) has a well-defined ex- 
pansion at k = 0, and R^(0,t) = 1 in agreement with 
the sum rule (1441). Thus, at zeroth-order we recover 



t 







Rl(x, t; qo) — ► Sjj(x — qo) and i?]?(fc, t) — ► 1. 

(66) 

This also applies to the large-scale limit k — > 0. On the 
other hand, Eq. lpS)) leads to the asymptotic behavior at 
large wavenumber K , 



K > 1 : ^(JsT) 



4^F 



Ai'(-wi) 



X COS 



V2 



K 



3/2 



^1 

V2 



AT 



(67) 



Thus, the propagator R^(K) decays somewhat more 
slowly than a Gaussian, with increasingly fast oscilla- 
tions, at high K. We display the propagator W'(K) in 
Fig. [21 We can note that the oscillations only appear in 
the very far tail where R^{K) is already negligible. 



10 



The density propagator R s of Eq. (|49|) directly follows 
from Rr using Eq. ({52"|) . In terms of dimensionless vari- 
ables, this gives 

R S (x, t; go) = t^tt R S (X - Qo), R s (k, t) = t R S (K), 



L(t) 



with 



(68) 



R S {X-Q ) = R*(X-Q ), R S {K) = R^(K). (69) 

Thus, the density propagator shows the same 
exponential-like decay at large wavenumbers or late 
times obtained in (|6"T)) for R^(k,t). 

For generic index n it is difficult to obtain the exact 
form of the velocity distribution p x (u), which is needed 
to derive the Eulerian propagators. However, as shown 
in [I?], HH , for — 1 < n < 1 we have the asymptotic large- 
velocity tails 



1 < n< 1, 



p(u) 



(70) 



where we do not specify prefactors and numerical factors 
in the exponents. This gives 



30 



00 



R^(x,t;q )~e- lx - qor+3/t2 . (71) 



As discussed in [3l| these behaviors are related to rare- 
event distributions. By contrast, the high-wavenumbcr 
tail in Fourier space depends on the details of the small- 
scale highly nonlinear processes and is more difficult to 
estimate. Note that the large-scale decay (|71[) depends on 
n, contrary to the Gaussian decay obtained in Eqs. (f55|) - 
(J5BJ) for the IR-class. 

As is clear from the expression ([BTj) and as discussed 
in section lTlIBl the Eulerian propagators R^(x, t; qo) and 
R s (x, t; qo) are governed by the one-point velocity distri- 
bution. We can note that various approximation meth- 
ods, based on resummation schemes or field-theoretic 
methods, have been recently devised to estimate such 
propagators in the case of the gravitational or Zel'dovich 
dynamics d, [TJ, [H, [l||, and they show a Gaussian or 
power-law decay that is set by the variance of the initial 
velocity, in agreement with the results obtained in sec- 
tion IIIIBl However, as pointed out in [ll|,|l9|], this decay 
does not express a true loss of memory associated with 
a relaxation towards some equilibrium but it is only due 
to the random advection of the density structures by the 
large-scale velocity effect (sweeping effect). For the case 
of white-noise initial velocity the initial velocity variance 
is not well defined but it becomes finite as soon as t > 
and the physics is the same. 

Thus, the strong decay of the Eulerian propagator is 
due to the random advection of the flow by the velocity 
field and it does not imply a true loss of memory for the 
structures themselves. For instance, shifting the linear 
density field Sl(x ) t) by a random uniform translation 
x 1— > x + aL(t), with a being distributed according to 
P(U) of Eq. fUTj) . would lead to the same decaying prop- 
agator even though structures in the density field remain 
unchanged. 



IV. LAGRANGIAN PROPAGATORS 
A. Definitions and linear regime 

In order to go beyond the apparent loss of memory due 
to the "sweeping effect" discussed in the previous sec- 
tion, one method is to work within a Lagrangian frame- 
work where the effect of uniform translations automati- 
cally vanishes. Thus, following [l7|, we consider the La- 
grangian propagator, R K (q, t; qo), associated with the La- 
grangian quantity n(q, t) defined by 



d dx 



whence n < 1, 



(72) 

where x(q, t) is the Lagrangian map that describes the 
trajectory of particle q. The upper bound, k < 1, is as- 
sociated with the fact that particles do not cross each 
other, so that x(q) is a monotonous increasing function. 
Thus, — k is the divergence of the Lagrangian displace- 
ment field x(3>*) — x (<?>*) — 3- It also describes the rel- 
ative expansion of infinitesimal mass elements, and from 
Eq. (f32|) it is related to the density field p(x,t) as 



p(x,t) = 



Po 



1 - K(q,t)' 



whence n(q, t) = 1 



Po 



P(q,ty 
. (73) 

with p(q,t) = p(x(q,t),t). Note that if x(q) is not 
monotonous (as would be the case for other systems such 
as the collisionless gravitational dynamics) one would 
need to keep the absolute value, \dq/dx\, for the Ja- 
cobian that appears in the expression of the density in 
([32| . which would violate the relationship (|73|) between 
k and p. In higher dimensions the relation ([75]) no longer 
applies since the density is determined by the determi- 
nant of the deformation matrix, \dxi/dqj\, whereas k is 
defined as the trace of this matrix. From a theoreti- 
cal point of view, it is more convenient to work with 
K(q,t) than with the Lagrangian density p(q,t) 7 defined 
as p{x(q, t),t), because the Lagrangian equations of mo- 
tion are usually most easily expressed in terms of the 
displacement field, x(3, t) = x(q,t) — q, whence in terms 
of its divergence — k (e.g. [l7l]). 

As is clear from Eqs. (TT2")) - (fT5|) , the variable k is not 
affected by uniform translations of the system, since it 
only depends on the divergence of the displacement field. 
Therefore, it does not suffer from the sweeping effect en- 
countered for the Eulerian variables in section ITTTI where 
we discussed Eulerian propagators. 

In the linear regime, the Lagrangian map is simply 
given by x = q + tuo(q), and we have from (|72p 



K L (q,t) = tK L o(q) with n La (q) = = ^ 



dq dq 2 



whence 



iiLo(k) = -fc 2 -0 Q (fc) = 6 L0 (k), 



(74) 



(75) 
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where we used Eq. (|51[) . Then, in a fashion similar to 
Eqs. (l4"2")l and (|49[) . we define the Lagrangian propagators 
R K (q,t;q ) and R K ^(q,t; q ) by 

R{q,t,q )-{ — — r), H ' v (q,t,q ) — { — — ). 

TJKLo(qo) Vipo{qo) 

(76) 

For systems that are homogeneous and isotropic these 
propagators only depend on the distance \q — qo\ and on 
time, and going to Fourier space we obtain 

pp. TDK Op. TDK 

= -k 2 R K , = "LfL. = < LP-. (77) 

oq z dqfi 

For systems in the IR-class, —3 < n < — 1, without IR 
cutoff, we take the limit where we are far from the refer- 
ence point to obtain an homogeneous system. 

In the linear regime, from Eqs. (jT4")) - (jT5"|) the Lagrangian 
propagators are given by 

Rl=t5 D (q-q ), R K L = t, (78) 



R^ =t8'^q-q Q ), R^ = -tk 2 . (79) 

The asymptotic behaviors apply to the early-time (t — > 
0) and large-scale (k — * 0) limits for the whole range 
— 3 < n < 1, as we shall check below. 



B. Relation with the shock mass function 

To compute the propagators beyond linear order, we 
first express n(q) in terms of the shocks built at a given 
time t, using the fact that all the matter at any time 
t > is located within shocks for the fractional Brown- 
ian motion initial conditions (|15p (34|. Therefore, any La- 
grangian point q belongs almost surely to a shock at time 
t > 0. Thus, for a given realization of the initial velocity 
field, let us note Xi the Eulerian position of the shock i, 
which gathers the particles coming from the Lagrangian 
interval [qi,qi + i[, with Xi < Xi+i and qi < q.i + i (we may 
choose to affect the index i — to the first shock to the 
right of x = 0). In the case —3 < n < —1, where shocks 
are dense, we only count shocks above a finite mass in- 
to obtain a discrete sum and we eventually take the limit 
m_ — * 0. Then, for qi < q < q^x we have x(q) = X{ 
whence n(q) = 1, whereas the jump from xi-\ to Xi 
at point qi gives a contribution — (x, — Xi-i)So{q — qi)- 
Therefore, we can write n(q) as 



K{q) = 1 - 



(Xi - Xi-i) S D (q - qi). 



(80) 



Next, in order to perform the functional derivative ([75]) 
we must evaluate the change of n(q) at linear order over 
a perturbation Sipo(qo). For smooth initial potentials 
ipo{q), the Lagrangian boundaries {qi^q^i} and the Eu- 
lerian location x% of the shock i are obtained from the 



geometrical construction (|30"f which writes as the four 
constraints 



Mli) = V^cMi) and ^ofe) = V'M^ ( 81 ) 



and 



V'ofe+i) = V Xi , Ci {q i+ x) and ^(ft+i) = V x Ci fe+i)- 

They express the condition that the first-contact 
parabola V XuCi simultaneously touches the curve at 
both points {qi,qi+i}, with a tangent slope. Applying 
the perturbation Sipo would give at point qi, 



t 



(Sqi - 8xi) + Sci, (83) 



H'oiqi) + ipo(qi)6qi = -(Sqi - Sxi 



(84) 



and similar relations at point qi+i, from which we can 
derive the changes 5qi,Sqi+i,Sxi, and Sci, in terms of 
Stpo. However, for the fractional Brownian motion initial 
conditions (|15p the changes Sqi and Sq^i vanish. 

Let us first consider the IR-class, —3 < n < — 1. Then, 
from the scalings (|20|) we can see that ip'o = ~ u o 1S finite 
but the second derivative ipQ = 9q is almost surely infinite 
(i.e. ijio has no regular second-order derivative; in the 
case n = —2 for instance, tp'oio) is a white noise). More 
precisely, from (|20|) we can see that the term ip'o(qi)Sqi in 
Eq. (f8"4"|) , that measures the change Aiip' of i(j' associated 
with a shift Sqi, is of order \Sqi\~^ n+1 ^ 2 , so that Sqi is 
of order \Sif>o\~ 2 ^ n+1 > (e.g., for n = — 2 where ip' (q) is 
a Brownian motion, Aiip' ~ vl&Zij and Sqi ~ |<5V>o| 2 )- 
Therefore, the changes Sqi and 5qi+\ are higher-order 
terms and do not contribute at linear order to Eq. (f83]) . 

On the other hand, for the UV-class, — 1 < n < 1, 
ipo(q) has almost surely no finite derivative and the 
change Aj^o of ipo associated with a s hift Sq.j, is of or- 
der | C 1 "")/ 2 (for instance, Aipo ~ vf^5i] for the case 
n = of white-noise initial velocity). Then, the tangent- 
slope constraints in Eqs.([8T j) -(f82 f do not apply and the 
changes Sqi and Sq^i are now non-perturbative, so that 
there is no term ip' (qi)8qi in Eq. (|53"l) . 

Therefore, for the whole range —3 < n < 1, the 
Lagrangian boundaries qi of the shocks do not change 
at linear order while the changes Sxi and Sci are ob- 
tained from Eq. (|83|) . and its companion at qi+\, setting 
Sqi = Sqi + i = 0. This gives 



Sxi = —t 



Sip (q i+1 ) - Sfojgj) 
q%+i - qi 



(85) 



Substituting into Eq. (|80|) yields 



Sn{q) = t ^2 s o(q-qi) 



#ofe+i) - Sipo(qi) 



Qi+i 



Sipo(qi) - 8tp {qi-i) 



qi-i 



(86) 
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and the definition (|T6[) gives 
R K ^(q,t;q') = t (£)M</ - ft) * 
& - gO - <fo [q% -q') 8 D (q i -q')-d D (q i -i-q') 



Qi+i - q% 



qi-i 



(87) 



This also reads as 
R K ^(q,t;q') = t (^M?"*) 

i 

(5 j (g - q' + A t q) + Sp jq-q'- A tg ) - 2fo (q-q') 
Xq 

where we noted A^q = qi+\ — qi and we used the homo- 
geneity and isotropy of the system (so that Aj_iq and 
Aiq have the same statistical properties). 

Let us define the shock mass function at time t, 
n(m, t)dm, as the mean number of shocks, per unit Eulc- 
rian or Lagrangian length (both functions are identical), 
with a mass in the range [m,m + dm]. Since the initial 
density is uniform and equal to po, the mass m of a shock 
is related to its Lagrangian size q (also called the "shock 
strength") bym = p q. Then, Eq. ([55|) writes as 



R K ^(q,t;q ) 



t 



dmn(m, t) 



8D{q-qo+m/p ) + 5 D (q-q -m/po) - 25 D (q-q ) 

x (»yj 

m/po 

We can note from this expression that this Lagrangian 
propagator satisfies the sum rule 



dqR K ^(q,t;q )= / dq R K ^{q, t; q ) = 0. (90) 



The propagator R K -^(q, t; q ) is singular at q = qo but for 
non-zero separation Eq. (159^1 simplifies as 



q^ qo 



R^(q,t;q )^tp o n(PQ } q - q °}' t \ 



q - 90 1 



(91) 



We can check that expression (|89|) only depends on the 
distance \q — qo I an d on time t, and it is clearly invari- 
ant through Galilean transformations. As expected, we 
can see that the Lagrangian propagator R K '^ provides 
a direct probe of the matter density field (which is di- 
rectly related to the structures of the displacement field, 
see Eqs. ([32|) and (|73|) ) as it can be expressed in terms 
of the shock mass function. Of course, the propagator 
R K ^ alone is not sufficient to fully characterize the den- 
sity field (which would require for instance all n-point 
correlation functions) but this represents a significant 
improvement over the Eulerian propagators described in 
section Hill In the linear regime, where the typical shock 
mass goes to zero while the shock mass function remains 
normalized to unity, 



/>oo 

I m . 

\ dm — n(m, t) = 1. 

Jo Po 



(92) 



or 
I 






Q-Qo 



I 




iog(|q QJ) 



FIG. 3: Upper panel: The Lagrangian propagator R K (Q — Qo) 
obtained for n = —2, in terms of dimensionless variables, 
from Eqs. (|95p and (|105p . The dashed lines are the asymptotic 
behaviors (|109[) - (|110[) . Lower panel: Same as upper panel but 
on a logarithmic scale. 



Eq. 
(179 



clearly converges towards the first expression 
Going to Fourier space we obtain 



R K ^(k,t) 



It I dmn(m,t) 
lo 



cos(km/pa) — 1 
m/po 



(93) 



which agrees with the second expression ((79)) in the limits 
t -> or k -> 0. Note that R K ^ and R K ^ are well 
defined over the whole range —3 < n < 1, thanks to the 
property (|92p which ensures that there is no divergence 
at low mass in Eq. (|93"l) . Thus, contrary to the Eulerian 
propagators studied in section HTT1 there is no qualitative 
difference between the UV-class and the IR-class for the 
Lagrangian propagator R K ^ (however, for the IR-class 
where the initial velocity field only shows homogeneous 
increments we must first take the limit of large distance 
from the reference point, to avoid boundary effects). 

The Lagrangian propagator R K follows from Eqs. ([77|) . 
(1591 and J93l) as 



R K (q,t- 1 qo) = t / dmn(m, t) 
Jo 

x \q-qo + m/p \ + \q-qo-m/p \ - 2\q-q \ 

2m/ p ' 1 ' 
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log(K) log(|Q QJ) 



FIG. 4: Upper panel: The Lagrangian propagator R K (K) in 
Fourier space, obtained for n = —2, in terms of the dimen- 
sionless wavenumber K, from Eqs. (|96|l and (I105|l , The dashed 
line is the asymptotic power-law behavior Lower panel: 

Same as upper panel but on a logarithmic scale. 



FIG. 5: Upper panel: The Lagrangian propagator R K (Q — Qo) 
obtained for n = 0, in terms of dimensionless variables, from 
Eqs. (|95p and (|112p . The dashed lines are the asymptotic cu- 
bic exponential behavior (|115[l . Lower panel: Same as upper 
panel but on a logarithmic scale. 



which also reads as 

R K (q,t-q ) = t [ dmn(m,t) (l-^—^i), 

Jp \q-qo\ V m /P0 J 

while we obtain in Fourier space 

#■(*,*) = 2* r^n{m,t) l - C ^ kmlp0 \ (96) 
Jo k-'m/po 

Note that Eqs.([M|)-([n5]) satisfy the constraint 
R K (q,t;qo) —> for \q — go I — * °°, hence no addi- 
tional integration constant appears when we use Eq. (|77|) 
to derive R K from R K ^ . 

As shown in [U, H3, 13 , the shock mass function shows 
the large-mass tail 

-3<n<l, m^>p L(t): n(m, t) ~ e -™" +3 /* 2 , (97) 

which is also related to the rare-event velocity tail ([70)) . 
Again we did not write numerical factors in the exponen- 
tial. At low masses, numerical simulations and heuristic 
arguments [3(1 HI] suggest the power-law tail 

-3<n<l, m<£p L(t): n(m,i) ~ t~ l m ( "" 1)/2 , 

(98) 



which has only been proved rigorously for the Brownian 
case n = — 2 [26|, [37], SH and for the white-noise case 
n = 0, HI- Then, Eqs.(H7J, (HJ), and ([77]), yield the 
large-separation asymptotic tails 

^^»1: R^(q,t-q n ) ~ R«(q,t;q ) ~ e -i™r +s /* a 

(99) 

which show the same behavior as the tails of the Eule- 
rian propagator (fTTj) . This is due to the fact that they 
are governed by the same rare events, associated with 
extreme fluctuations in the initial velocity field [3l|. At 
small separations, Eqs. (f9"8"|) and (f9"Tj) give 

\q-q \<L(t) : R«>*{g, t; q ) ~ \q - « | ( "- 3)/2 (100) 
R K (q,t;q ) [-i? K (0, t; 0)] ~ \q - q Q \ {n+1)/2 , (101) 

with, for — 1 < n < 1, 

/» CO 

R K (0,t;0)=t dmn(m,t) cx t ( ™ +1)/( " +3) . (102) 
Jo 

The term i? K (0, i; 0) is only present in Eq. (|10ip for -1 < 
n < 1 (for —3 < n < — 1 the propagator R K (q,t;qo) 
diverges for q — » q a as |g — g |^™ +1 ^ 2 )- Note that the 
time-dependence (|102p obtained for — 1 < n < 1 is such 
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1 o 
K 




log(K) 



FIG. 6: Upper panel: The Lagrangian propagator R K (K) in 
Fourier space, obtained for n = 0, in terms of the dimension- 
less wavenumber K, from Eqs. (|96[l and (|112|l . The dashed 
line is the asymptotic power-law behavior (|116[) . Lower panel: 
Same as upper panel but on a logarithmic scale. 



that we recover the linear behavior ([75)1 on large scales, 
as we have R K (0,t;0)L(t) <x t. For — 3 < n < — 1 the 
linear time-dependence (|78p is also recovered on large 
scales from Eq.flTUI]) through L{t) { - n+1 ^ 2+1 oc t. On the 
other hand, for all —3 < n < 1, the time dependence 
vanishes in the leading-order terms in the right hand side 
of Eqs.([100]) and (flUT]) . 

In Fourier space, from Eqs. (|9"3"l) , (|96p . we recover at low 
k (k ->• 0) the asymptotic limits (J7HJ> and 1(75]). At high 
A:, the small-separation singularities (|100|) - (|101[) give rise 
to the time-independent power-law tails 



k^>L{t)- 1 : R K ^(k,t) ~ Z^ 1 -™)/ 2 , (103) 
R K (k,t) ~ fc-("+ 3 )/ 2 . (104) 



Thus, R K '^(k,t) grows at high k whereas R K (k,t) slowly 
decreases, in a time-independent fashion, over a range of 
wavenumbers that grows with time since L{t)~ l goes to 
zero at late times. 

Note that, contrary to the Eulerian propagators dis- 
cussed in section IIIIB1 the Lagrangian propagators are 
well defined over the whole range — 3 < n < 1. As ex- 
plained in section IIV A[ the Lagrangian quantity n(q, t) 
is not modified by uniform translations of the system, so 
that the propagators R K '^ and R K are not sensitive to 



long- wavelength modes of the initial velocity and remain 
finite for — 3 < n < —1 as we push a possible IR cutoff 
to infinity. 



C. Brownian and white-noise initial velocity 

For the Brownian case, n = —2, the shock mass func- 
tion is explicitly known [2ol |37| , 



n = -2 : N(M) = 



Af-3/2 



— M 



(105) 



where we introduced the dimcnsionlcss scaling variables 
m , 1 



M 



n(m, t) 



Po L(ty " v "" v Po L(t)2 

Introducing the dimensionless propagators, 
t 



N(M). (106) 



R K (q, t; q Q ) = j^R^Q - Qo), R K (k, t) = t R R (K), 

(107) 

we obtain from Eq. (|9~lj) 



&R K r^ _ N(Q) 
12 w> ~ n 



-5/2 



,-Q 



dg 2 



whence 



Q<1 : R K (Q) 
Q>1 : R K {Q) 
as well as R K (0) = 1 and 



3V^Q' 



(108) 

(109) 
(110) 

(111) 



For the white-noise case, n = 0, the shock mass func- 
tion is again explicitly known [U El[ , 

f+ icc ds' e~ s ' M 
n = 0: N(M) = 2M / _ ___ 
' 27n Ai(s') 2 



>/ — 200 

+ 200 



ds 



ZTTl Al(s) 



using the dimensionless scaling variables Q106P . It shows 
the asymptotic tails [4l| 

M»l : N{M) ~ 2^M b/2 e- bJlM - M3/12 , (113) 
M<1 : N(M) ==, (114) 



which agree with Eqs. (f97|) - ([98|) . In terms of the dimen- 
sionless propagators (|107j) we obtain R K (0) ~ 1.674 and 
the asymptotic tail 



Q»l: R K (Q) ~ 32v / 7rQ -5 ' 2 e 



5/2 -W!Q-Q 3 /12 



(115) 
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n — O 



i = — 2 




FIG. 7: Upper panel: The Lagrangian propagator R K (T) 
as a function of the dimensionless time T = t/tk, at fixed 
wavenumber k. The dot-dashed line corresponds to Brown- 
ian initial velocity (n = —2) and the solid line to white-noise 
initial velocity (n = 0). Lower panel: Same as upper panel 
but on a logarithmic scale. 



as well as R K (0) = 1 and 



(116) 



We show in Figs. [3][6] the propagators — Qo) and 

R K (K) obtained for the two cases n = —2 and n = 0. 



D. Evolution with time 

As shown in the previous section, for a fixed wavenum- 
ber k we recover at early times the linear growth ([78]) - 
([79]). whereas at late times the Lagrangian propagators 
saturate into the nonlinear regime, as seen in Eqs. (|103p - 
(|104| . To emphasize the evolution with time, we define 
the characteristic time t kl associated with wavenumber 
k, which measures the time when scales ~ 1/k enter the 
nonlinear regime (i.e. K = 1), by 



kL(t k 



1, whence t k oc k 



-(n+3)/2 



(117) 



where Lit) was introduced in Eq. (|39p . and we introduce 
the associated dimensionless time-variable T, 

t 

tk 



Then, we write the propagator R K (k,t) of Eq. (|107p as 

R K (k,t) = t k R K (T) whence R K (K) = ^R K (T). (119) 

As noticed above, the propagator R K (T), which describes 
the time-dependence of R K (k,t) at fixed k, shows the 
early- and late-time asymptotics 



T« 1 : R K (T) ~ T, T 



T 



whence T = K (rl+3)/2 . 



(118) 



a : R K (T) -> constant. 

(120) 

We obtain i? K (oo) = 4^/3 for n = -2 and R K (oo) = 
2y2 for n = 0. We display in Fig. [TJthe scaling propaga- 
tor R K (T) obtained for n = —2 and n = 0. 

The behaviors (|120[) explicitly show that the La- 
grangian propagators first grow linearly with time un- 
til they reach the nonlinear regime where they saturate, 
by contrast with the Eulerian propagators, which show 
exponential-type decays as seen in section Hill The fact 
that the former show no decay at late times implies that 
there is no real loss of memory. 



V. CONCLUSION 

In this article we have studied the Eulerian and La- 
grangian propagators which are obtained within the ad- 
hesion model (i.e. Burgers dynamics), focusing on the 
one-dimensional case with power-law initial conditions 
that develop a self-similar evolution. We have derived 
some general and asymptotic results for linear power- 
spectrum index n in the range — 3 < n < 1, as well 
as complete explicit expressions for the two representa- 
tive cases n = — 2 (Brownian initial velocity) and n = 
(white-noise initial velocity). In particular, we note that 
the range — 3 < n < 1 can be split over an "IR class", 
—3 < n < — 1 , governed by long wavelengths, and a "UV 
class", — 1 < n < 1, governed by small wavelengths. 

We find that Eulerian propagators can be expressed in 
terms of the one-point Eulerian velocity probability dis- 
tribution. This clearly shows that they are sensitive to 
long-wavelength modes of the velocity field, which leads 
to a "sweeping effect" as small-scale structures can be 
moved over large distances without significant distortions 
if most of the power is stored at very low wavenumbers. 
This yields a strong exponential-like decay at high-fc of 
Eulerian propagators, which does not imply a strong loss 
of memory of the system as it is due to this random ad- 
vection. In the IR-class, the Eulerian propagators show 
an universal Gaussian decay (i.e. independently of n) of 

the form e "o' in Fourier space. However, if there 
is no infrared cutoff the factor a^ diverges and the Eule- 
rian propagators vanish as soon as t > 0. In the UV-class, 
Eulerian propagators strongly depend on the index n. In 
particular, for n — we obtain in Fourier space an oscil- 
latory decay, with an amplitude of the form e~ tk 1 . In 
both IR and UV classes, being related to the one-point 
Eulerian velocity probability distribution Eulerian prop- 
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agators are not sensitive probes of the structures of the 
underlying density field. 

For more complex dynamics, such as cosmological 
gravitational clustering, Eulerian propagators can no 
longer be written in terms of the velocity distribution 
in such a direct manner, but they remain governed by 
the same sweeping effect and resummation schemes typ- 
ically show a decay at high wavenumbers and late times 
[1, EH, E3- Even when this effect is modest (i.e. 
there is not much power at low k) Eulerian propagators 
are unlikely to provide good probes of the density field, 
its relation to the one-point velocity distribution being 
rather loose. 

Next, we have shown that Lagrangian propagators can 
be expressed in terms of the shock mass function (which 
corresponds to the halo mass function in the cosmo- 
logical context) as shown explicitly in Eqs. (|89|) . (f93|) 
and (|96[) . Therefore, Lagrangian propagators are much 
more closely related to the properties of the density field. 
Moreover, they show the same properties for both the IR 
and UV classes, as could be expected from the Galilean 
invariance of the equations of motion. Whereas on large 
scales and low k they grow with time, in agreement with 
the linear regime, on small scales and high k they sat- 
urate to constant values, with a power-law dependence 
on scale or wavenumber that is set by the initial index 
n. This strong memory of the system marks a sharp 
contrast with the decay of different-time Eulerian statis- 
tics. In higher dimensions we still expect the Lagrangian 
propagators to show power-law tails over wavenumber k 



in the highly nonlinear regime, with a slope that is again 
related to the low-mass tail of the shock mass function, 
but there can also be a power-law decay over time. We 
leave such a study to future work. 

Even though in more complex dynamics Lagrangian 
propagators are unlikely to be expressed in terms of the 
density field in such a direct manner (i.e. through the 
mass function of bound objects), they should provide a 
sensitive probe of the behavior of the density field and of 
the relaxation processes at work. As seen in this article, 
a great interest of the adhesion model is to provide exact 
results in a nontrivial case which shares some key prop- 
erties with cosmological gravitational clustering. This 
allows us to obtain some insight in the behavior of com- 
plex quantities, such as propagators or response func- 
tions, and the processes which they probe. Then, this 
can serve as a guide to decipher the processes associated 
with more complex systems. On a quantitative level, it 
remains to be seen whether this can be used to build 
for instance efficient ansatze (in the manner of the halo 
model [i| or the stable-clustering ansatz [l|) that would 
help building an accurate description of the system. 
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